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Abstract. Precision measurements of the cosmic microwave background by WMAP are believed to have estab- 
lished a flat A-dominated universe, seeded by nearly scale-invariant adiabatic primordial fluctuations. However 
by relaxing the hypothesis that the fluctuation spectrum can be described by a single power law, we demonstrate 
that an Einstein-de Sitter universe with zero cosmological constant can flt the data as well as the best concordance 
model. Moreover unlike a A-dominated universe, such an universe has no strong integrated Sachs- Wolfe effect, so 
is in better agreement with the low quadrupole seen by WMAP. The main problem is that the Hubble constant 
is required to be rather low: Ho — 46 km/s/Mpc; we discuss whether this can be consistent with observations. 
Furthermore for universes consisting only of baryons and cold dark matter, the amplitude of matter fluctuations 
on cluster scales is too high, a problem which seems generic. However, an additional small contribution {Qx ~ 0.1) 
of matter which does not cluster on small scales, e.g. relic neutrinos with mass of order eV or a 'quintessence' with 
U) ~ 0, can alleviate this problem. Such models provide a satisfying description of the power spectrum derived 
from the 2dF galaxy redshift survey and from observations of the Ly-a forest. We conclude that Einstein-de Sitter 
models can indeed accommodate all data on the large scale structure of the Universe, hence the Hubble diagram 
of distant Type la supernovae remains the only direct evidence for a non-zero cosmological constant. 
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1. Introduction 

Measurements of cosmological parameters with reasonable 
accuracy are essential both to establish a robust picture of 
the standard Big Bang cosmology, and to provide insights 
into the fundamental processes, far beyond the Standard 
Model of particle physics, which determined its initial con- 
ditions. Since the pioneering work of Hubble, it has been 
recognised that cosmological tests based on astrophysical 
arguments can suffer from large systematic biases. Clearly 
one should as far as possible use methods which do not 
depend explicitly on assumptions concerning complex as- 
trophysical phenomena. 

In this respect, measurements of anisotopies in the 
Cosmic Microwave Background (CMB) appear to offer the 
most promise for accurate determination of cosmological 
parameters, thanks to the high control possible on system- 
atic errors. Since the epochal discovery of primordial fluc- 
tuations on large angular scales by COBE (Smoot et al. 
1992), this field has witnessed a renaissance. The first de- 
tections of fluctuations on degree scales (Netterfield et al. 
1995; Scott et al. 1996) provided tantalizing evidence for 



Send offprint requests to: Alain. BlcOichardSast . obs-mip . fr 



the flatness of the Universe (e.g. Lineweaver et al. 1997). 
The unambiguous detection of the first and second acous- 
tic peaks in the angular power spectrum (de Bernardis et 
al. 2000; Hanany et al. 2000; Halverson et al. 2002) has 
confirmed this result. Taken together with studies of large- 
scale structure (LSS) in the universe, these observations 
have also confirmed the overall picture of structure forma- 
tion through gravitational instability. The recent results 
obtained by WMAP represent a further major advance in 
the field. For the first time, measurements of cosmological 
parameters are being quoted with uncertainties of a few 
per cent, opening up the anticipated era of 'precision cos- 
mology'. Our intention here is to examine whether such 
determinations are in fact robust or depend crucially on 
underlying assumptions or 'priors'. Specifically we wish to 
test whether a cosmological constant. A, is really required 
by observations of the CMB and LSS, independently of the 
indications from the SN la Hubble diagram. We will do 
so by confronting Einstein-de Sitter (E-deS) models with 
the same observations. It turns out that with a different 
assumption concerning the spectrum of primordial fluctu- 
ations generated by inflation, such models can fit the data 
even better than models with non-zero A. 
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2. What do the Q measurements imply? 

The physics of passive linear perturbations in the early 
Universe is well understood and therefore their evolution 
can be computed accurately (see Hu & Dodelson 2002). 
This is the basic reason why precise measurements of the 
angular structure of the imprints left on the CMB by 
primordial fluctuations can provide accurate information 
on cosmological parameters. The ingredients necessary to 
compute the amplitudes of the multipole moments (C^) 
specifying the angular power spectrum are both the na- 
ture and spectrum of the primordial fluctuations (presum- 
ably arising from inflation), as well as specification of the 
various contents of the universe which contribute to its 
density and pressure. The imprint of a specific parameter 
can be direct, through the infiuence on the dynamics of 
acoustic oscillations before the epoch of last scattering (as 
for the baryonic content for instance), as well as indi- 
rect, through the effect on the angular distance to the last 
scattering surface. A non-zero A affects the Ces primar- 
ily through the distance effect (Blanchard 1984). There 
are also more subtle effects, such as the integrated Sachs- 
Wolfe (ISW) effect which contributes at a much weaker 
level. However such effects are harder to identify, as they 
can easily be mimicked by a non-trivial primordial fiuctu- 
ation spectrum. 

The first studies of the generation of density pertur- 
bations during inflation established (see Linde 1990) that 
for the simplest models involving a single 'inflaton' field, 
the spectrum is close to the Harrison-Zel'dovich (H-Z) 
scale-invariant form, P{k) oc k" with n = 1, which had 
been proposed earlier on grounds of simplicity. Thus the 
H-Z spectrum became a standard input for calculations 
of CMB anisotropics and the growth of LSS, e.g. in the 
standard cold dark matter (SCDM) model (Davis et al. 
1985). In fact there are significant corrections to a H- 
Z spectrum even in single-field models, in particular the 
spectrum steepens logarithmically with increasing k (de- 
creasing scale) as the end of inflation is approached. This 
is usually accomodated by considering a 'tilted' spectrum 
with n < 1, although it should be noted that the index 
n is scale- dependent for any polynomial potential for the 
inflaton, and is constant only for an exponential potential 
(power-law inflation). Moreover n can be close to, and even 
exceed, unity if inflation ends not through the steepening 
of the inflaton potential but, for example, due to the dy- 
namics of a second scalar field (hybrid inflation). In such 
multi-field models, the spectrum may not even be scale- 
free since features can be imprinted onto the spectrum, e.g. 
when the slow-roll evolution of the inflaton is interrupted 
by other background fields undergoing symmetry-breaking 
phase transitions (Adams, Ross & Sarkar 1997b). 

The expectations for the spectral index n{k) in various 
inflationary models has been reviewed by Lyth & Riotto 
(1999). Even small departures from scale-invariance can 
be quite significant for LSS formation. For example after 
the SCDM model was foimd by Efstathiou et al. (1992) 
to be in conflict with the observed power spectrum of 



galaxy clustering (in being unable to simultaneously re- 
produce the abundance of rich clusters (quantifled by the 
variance erg in a top-hat sphere of radius 8h~^ Mpc;) and 
the COBE measurement of fluctuations on the scale of 
Hg^ ~ 3000/1" 1 Mpc), it was noted by White et al. 
(1995) that invoking a tilted spectrum with n ~ 0.9 could 
save the model. Interestingly enough, such a spectrum 
arises from natural supcrgravity inflation,^ where the lead- 
ing term in the potential is cubic in the field (Ross & 
Sarkar 1996; Adams, Ross & Sarkar 1997a). This yields 
n = {N-2)/{N + 2), where < 57-Hln {k~^/H^^) is the 
number of e-folds of expansion from the end of inflation, 
taking the inflationary energy scale to be ^ 10^^ GeV as 
required by the normalization to COBE, and the reheat 
temperature to be < 10^ GeV to avoid the thermal grav- 
itino problem (see Sarkar 1996). If the inflationary scale is 
signiflcantly lower (German, Ross & Sarkar 2001) and/or 
if there was a late epoch of thermal inflation (Lyth & 
Stewart 1996), then our present Hubble radius may have 
exited the horizon only ~ 20 — 30 e-folds from the end of 
inflation. This yields a spectral index as low as n ~ 0.8 on 
cosmologically observable scales since n ~ 1 — A/N in this 
model. 

It is thus clear that the primordial spectrum may not 
have a trivial form and lacking a 'standard model' of in- 
flation, it is necessary to consider a wide range of pos- 
sibilities. Furthermore such complex spectra could po- 
tentially confuse cosmological parameter estimation from 
CMB data (e.g. Kinney 2001). This was explicitly demon- 
strated by Barriga et al. (2001) using the COBE and 
BOOMERanG data (de Bernardis et al. 2000) for the case 
of a primordial spectrum with a step-like feature at a scale 
k 0.1 h/Mpc, as in double inflation (Silk & Turner 
1987). In this paper we investigate the flexibility in the 
determination of cosmological parameters using the much 
more precise WMAP data, when the usual hypothesis of a 
single power law spectrum is relaxed. Since the power-law 
index is related to the slope and curvature of the inflaton 
potential, it can change suddenly e.g. if the mass of the in- 
flaton changes through its coupling to a background field 
which undergoes spontaneous symmetry breaking during 
inflation (Adams et al. 1997b) 

The possible detection of a non-zero A through mea- 
surements of the Hubble diagram of distant Type la super- 
novae (Riess et al. 1998; Perlmutter et al. 1999) is among 
the most significant developments in modern cosmology 
(see Peebles and Ratra, 2003), and has led to the estab- 
lishment of the 'concordance model' with ^ 0.7 and 
ilm ^ 0.3 (Bahcall et al. 1999). This is quite consistent 
with the WMAP data for an assumed power-law primor- 
dial spectrum; for a spatially flat universe the cosmologi- 
cal parameters are determined to be: 51,„/i^ = 0.14 ±0.02, 
n^h^ = 0.024 ± 0.001, n^h^ < 0.0076, h = 0.72 ± 0.05, 
n = 0.99 ± 0.04 and r = 0.166tg;g^f (Spergel et al. 2003). 

^ Technically 'natural' means that the flatness of the poten- 
tial is protected by a symmetry — here the shift symmetry of 
a Nambu-Goldstone mode (Freese, Freiman & Olinto 1990). 
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This agreement has led to to the widespread behcf that 
the ACDM concordance model is now established to high 
accuracy through CMB measurements. However we wish 
to illustrate that this agreement is crucially dependent 
on the underlying assumptions concerning the primordial 
power spectrum and that CMB data do not yet indepen- 
dently require a non-zero A. Before addressing this issue, 
let us first assess the status of the concordance model. 



variance (AT| = £{i + l)Ce/27r) of AT| = 154 ± 70 iiK"^ 
(Bennett et al. 2003), so one might conclude that e.g. an 
expected AT2 — 350 ^K^ is discrepant by about 2a, which 
would be at the 95% c.l. for a gaussian distribution. In 
properly evaluating this probability however one should 
take into account the foreground removal technique and its 
uncertainty, in order to determine the likelihood distribu- 
tion. Such a distribution is likely to be non-Gaussian and 
it is therefore possible that the estimation of the goodness- 
of-fit for the concordance model might be significantly im- 
proved in this outlying region. If one instead considers the 
best-fit concordance model derived by Spergel et al. (2003) 
to be the true description of the CMB sky, the probability 
of observing a low quadrupole can be directly estimated. 
The log-likelihood of Cf^°° can be well approximated by 



(Bartlett et al. 2000): 



21n/:(C)*^™) =5x/,kyX 
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3. Is the concordance model actually concordant? 

3.1. Concordance with WMAP data 

As mentioned already, SCDM with ilm = 1, h — 0.5 and 
n — 1 was found to disagree with the shape of the APM 
galaxy correlation function (Maddox et al. 1990), as well 
as the high baryon fraction measured in clusters, which 
together with the baryon fraction inferred from primor- 
dial nucleosynthesis arguments, implied a lower matter 
density of fim ~ 0.3 in agreement with local dynamical 
estimates (White et al. 1993). The subsequent measure- 
ments of CMB fiuctuations on degree scales however re- 
quired the universe to be spatially flat and ruled out such 
a low density matter-dominated universe (Lineweaver & 
Barbosa 1998). Thus the possible detection of cosmic ac- 
celeration in the Hubble expansion of distant supernovae, 
implying a cosmological constant with JIa ~ 0.7, was ea- 
gerly seized on as a mean of reconciling the CMB and 
LSS data. However although this concordance model is 
consistent with most cosmological observations, its first 
precision test has come with the WMAP data. The agree- 
ment of the concordance model with data as summarised 
by the WMAP team appears impressive (Spergel et al. 
2003). However, there are two facts to keep in mind. First 
the global on the temperature autocorrelation (TT) 
power spectrum is rather poor — the probability that the 
model fits the data is only 3%, so strictly speaking the 
model is rejected at the ~ 2a confidence level! However 
given possible remaining systematics eflfects not yet ac- 
counted for, the WMAP team concluded that this should 
not be considered as a serious problem for the concordance 
model. Allowing for 'running' of the spectral index with 
scale improves the fit somewhat; the data suggest that 
n > 1 on the largest scales and n < 1 on small scales. 

There is another aspect of the WMAP data that is even 
more puzzling, viz. the amplitudes of the low C^s, par- 
ticularly the quadrupole, is rather small compared to the 
expectation in the concordance model where the large cos- 
mological constant should boost the anisotropy on large 
angles. It is well known that the cosmic variance is high 
on such large angular scales and that Galactic foreground 
subtraction introduces further uncertainties. Spergel et al. 
(2003) concluded, from Monte Carlo realisations following 
two different methods, that the low signal on large scales 
cannot be obtained in over 99% of the cases. 

However, statistical inferences from the quadrupole 
amplitude have to be handled with caution. The measured 
value of Qrms — \/ (5/47r)C2 = 8 ± 2 /iK corresponds to a ^ |http : //lambda . gsf c. nasa . gov. 



where the noise is K (= 3.4/j,K^ for ^ = 2 as quoted in the 
WMAP data release^), the sky coverage is /sky = 85%, 
Ci is the measured amplitude and C'f^°° is the value of 
the best-fit ACDM model (1204 for t = 2). Given 
this approximation, one can retrieve the probability dis- 
tribution of C2 and thus estimate the chance of observ- 
ing a low value, following Douspis et al. (2003). For 
Qrms — 8/iK this is 4.6%, in agreement with the estimate 
of Tegmark, de Oliveira-Costa & Hamilton (2003) by a 
different method (varying the cut CMB sky and consider- 
ing the best running spectral index model). This indicates 
that the quadrupole is an outlier at most at the 2a level 
(see also Gaztaiiaga et al. 2003, Efstathiou 2003b). 

3.2. Concordance with astronomical data 

As discussed already, the concordance model has been 
built up over time in order to match observations, thus 
its a posteriori agreement with much of the LSS data is 
not a test. Of course as the quality of data improves the 
model will be further tested, although the number of free 
parameters provides some room for adjustment. 

Interestingly enough, WMAP has thrown new light on 
the masses of galaxy clusters and therefore on the inferred 
baryon fraction which, it had been argued (White et al. 
1993), indicates a low matter density universe. There has 
been some controversy in recent years concerning the ac- 
tual masses of X-ray emitting clusters, which are deter- 
mined by two different methods. One is the application 
of hydrostatic equilibrium, while the second uses mass- 
temperature relationships derived from numerical simula- 
tions. Systematic differences between the two methods are 
significant (Markevitch 1998; Roussel, Sadat & Blanchard 
2000) and this translates into an appreciable difference in 
the derived value of a^ (Reiprich & Bohringer 2002; Seljak 
2002). Furthermore similar differences arise from the use of 
different theoretical mass functions. The Sheth & Tormen 
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(1999) expression is recognised as providing a satisfactory 
fit to the mass function obtained from numerical simu- 
lations. Using this for an Q^-^ = 0.3 universe, one finds 
CTg — 0.86 corresponding to the high mass estimates from 
numerical simulations, and erg = 0.68 for the low mass 
estimates from hydrostatic equilibrium (Vauclair et al. 
2003). Clearly the WMAP measurement of cts = 0.9 ±0.1 
(Spergel et al. 2003) favours the high mass estimates. The 
implied baryon fraction (including stars) in clusters is then 
rather low, slightly below 9% (for h = 0.7). This conflicts 
with the universal baryon fraction of 15% required for the 
best concordance model fit to the acoustic peaks! We note 
that this discrepancy disappears for as ^ 0.7, which is 
about a 2a deviation from the WMAP determination. 

4. Acceptable Einstein-de Sitter models 

Let us now examine whether it is possible to obtain an 
acceptable CMB power spectrum in an E-deS universe. 
Clearly to do this we have to deviate from the assump- 
tions of the ACDM concordance model regarding the pri- 
mordial power spectrum. Indeed the WMAP team (Peiris 
et al. 2003) have already noted that the model fit can be 
significantly improved (particularly to the outliers at £ = 
22, 40 and 200) by allowing for oscillations in the primor- 
dial spectrum (Adams, Cresswell & Easthcr 2001) such as 
might be induced by phase transitions occuring during in- 
flation (Adams, Ross & Sarkar 1987b). However if we are 
not to introduce too many new parameters, the simplest 
modification that can be introduced is perhaps to con- 
sider a change in the slope of the spectrum at a particular 
scale. It is important in this respect to notice that the first 
and second acoustic peaks span a rather limited range of 
scales, £ ~ 150 — 600, while the rising part of the first 
peak covers a much bigger range, £ ^ 2 — 200. Without 
advocating any specific scenario, it is interesting to exam- 
ine how a model with different power law indices for the 
primordial fluctuations in these two regions compares to 
the observational data. We therefore focus on models with 

p.,^ ^ / Aifc«i forfc < fcl, 
^ ' \ Aafc"^ for k > fcl, 

with a continuity condition (Aifc"^ = A2fc"^). We cal- 
culate the CMB power spectrum using the CAMB code 
(Lewis et al. 2000) and use the WMAP likelihood code 
(Verde et al. 2003) to determine the quality of the fit. 

As we are primarily interested in examining the pos- 
sible constraints on the cosmological constant, we have 
restricted our search to models with fl^ = but allow 
a reasonable range for other cosmological parameters (in- 
cluding the optical depth t to last scattering). The best 
model we find has h = 0.46, cjb = ^hh"^ = 0.019, t = 0.16, 
fcl = 0.0096 Mpc-i, ni = 1.015, riz = 0.806. As seen in 
Fig. n the calculated power spectrum does very well in 
fitting the WMAP data and other observations at high £. 
Interestingly enough a preferred scale of fc ~ 0.01 Mpc~^ 
was also found by Mukherjee & Wang (2003) in attempt- 
ing to reconstruct the primordial spectrum in the context 



of a ACDM model (although Bridle et al. (2003) did not 
detect this using a different method). We wish to empha- 
sise that inspite of having one additional parameter, our 
model has a better (on the scalar Ci) than the best 
concordance model, because of its lower amplitude at low 
I (the of the TE spectrum being identical). In particu- 
lar, the mean quadrupole C2 has an amplitude of 844 /iK^, 
which has a 13% probability of yielding Qrms = 8 /iK. The 
reason of course is that E-deS models do not produce ISW 
effects as high as in flat models with low matter density 
and a large cosmological constant. This is arguably a sim- 
pler way to accomodate the observed low signal at Z < 20, 
than to invoke new physics (e.g. Spergel et al. 2003; Uzan, 
Kirchncr & EUis 2003a; Efstathiou 2003a, Contaldi et al. 
2003, Chne, Crotty & Lesgourgues 2003). 
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Fig. 1. The temperature power spectrum for the 
best-fit power-law ACDM model (dotted black line) 
from Spergel et al. (2003), and for our broken-power- 
law model with JIa = (solid blue line), com- 
pared to data from WMAP and other experiments 
at small scales — CBI f Pearson e t al. 200 2|), AC BAR 
IIKuo et al. 2002)l, BOOM ERanG fRuhl et al. 2002|l and 
VSA ( |Grainge et al. 2003| ). Note the Hnear scale in I for 
I > 200. 



4.1. Why this E-deS model is in difficulty 

It might appear that our E-deS model with = 1, 
r^A = and h = 0.46 must be in conflict with a num- 
ber of astronomical observations. However, several of these 
observations which in fact support the ACDM concor- 
dance model have been questioned. For example, the ob- 
served mass-to-light function from galaxies to superclus- 
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ters yields 17^ = 0.16 ±0.05 (Bahcall et al. 2000) but such 
observations are mostly local and possibly untypical. The 
abundance of high redshift X-ray selected clusters, a global 
test rather than a local one, systematically leads to high 
values of fim, well above the best WMAP value (Henry 
1997; Sadat, Blanchard & Oukbh 1998; Viana & Liddlc 
1999; Borgani et al. 1999; Reichart et al. 1999; Blanchard 
et al. 2000). This points in the same direction as the ob- 
served absence of any large-scale correlations between the 
COBE map of the CMB and the HEAO-1 map of the 2-8 
keV X-ray background which provides an interesting up- 
per limit of < 0.6 at 95% c.l. (Boughn, Crittenden & 
Koehrscn 2002). ^ 

The only direct evidence so far for a cosmological con- 
stant comes from the Hubble diagram of distant Type la 
supernovae, a method which relies on the standard candle 
hypothesis and on empirical corrections to the observed 
peak magnitudes on the basis of the observed decay times. 
Such corrections are essential for reducing the scatter in 
the data sufficiently so as to allow significant cosmologi- 
cal deductions. However there are systematic differences 
in the corrections made for the same objects by the two 
groups (Leibundgut 2000) which raises legitimate concerns 
about their validity. Moreover, the distant SN la appear to 
be significantly bluer than the nearby sample, suggesting 
that the derived reddening may have been underestimated 
(Leibundgut 2001). Rowan-Robinson (2002) has argued 
that when extinction and the luminosity-decay time rela- 
tion are treated in a self-consistent way, the significance 
of the evidence for positive A is much reduced. 

A Hubble constant of Hq — 46 km/s/Mpc would seem 
to be completely inconsistent with the Hubble Key Project 
(HKP) determination of 72±8 km/s/Mpc (Frcedman et al. 
2001). However there are some details of this work which 
might merit reexamination; 

1. Rowan- Robinson (2003) finds that if a more sophisti- 
cated local flow model is used than that of Mould et 
al. (2000), there is a reduction of about 2 km/s/Mpc 
in Hq. 

2. The method of combining the data used by Freedman 
et al. (2001), viz. estimating Ho by different methods 
and then combining the results, docs make the out- 
come somewhat vulnerable to Malmquist bias. Rowan- 
Robinson (2003) finds that this can also result in over- 
estimation of Hq by about 2 km/s/Mpc. 

^ Very recently, Boughn & Crittenden (2003) claim to have 
detected a large-scale correlation between the WMAP and 
HEAO-1 maps at the 2.4 — 2.8a level, which is surprising since 
the WMAP and COBE maps are quite consistent. They also 
find a correlation between WMAP and the NVSS survey of 
radio galaxies at the 1.8 — 2.3a level. Nolta et al. (2003) con- 
firm the latter finding and reject a A = universe at 95% c.l. 
However Myers et al. (2003) have found significant contaminar- 
tion of the WMAP data by NVSS sources which can account 
for the observed correlation. Fosalba & Gaztanaga (2003) also 
find a correlation, twice as strong as the expected signal, be- 
tween WMAP and the APM galaxy survey. 
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k = 0, z > 0.05 




luminosity distance (Mpc) 

Fig. 2. Velocity versus luminosity-distance for Type la 
supernovae (filled circles), S-Z clusters (open circles) and 
gravitational lens time-delay systems (filled triangles), 
with z > 0.05. All curves shown correspond to flat models 
and are labelled with the Hubble parameter in km/s/Mpc. 

3. Between the summary paper of Mould ct al. (2000) 
and that of Freedman et al. (2001), a change in the as- 
sumed I-band period-luminosity relation resulted in an 
increase of the extinction values in the Cepheid pro- 
gram galaxies, and hence in an average reduction in 
the distance scale of 4%. This was almost exactly can- 
celled out by correction for the effects of metallicity 
on Cepheid distances, of about the same magnitude 
in the opposite direction. However it is possible that 
these effects do not quite cancel, e.g. if the extinction 
in the program galaxies has been overestimated and/or 
the correction for the effects of metallicity have been 
underestimated. It would be highly desirable to extend 
observations of the Key Project Cepheids into the in- 
frared to assess these effects. 

4. The assumed extinction in the LMC may be slightly 
on the low side compared with estimates for hot stars 
in the LMC by Zaritsky (1999). 

5. Finally, it would be highly desirable to confirm the 
assumed distance to the LMC, preferably by direct ge- 
ometric methods. 

All these effects are likely to be small, but it is possible 
that they may combine in the same direction to signifi- 
cantly reduce Hq. 

It is interesting that methods which are largely in- 
dependent of the LMC and Cepheid distance scales, do 
tend to give significantly lower values for Hq. For instance, 
Sunyaev-Zeldovitch (S-Z) distances to 41 clusters give a 
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value of 54^3 km/s/Mpc in an E-deS universe (Reese et 
al. 2002); furthermore any clumping of the X-ray emit- 
ting gas would lower the actual value by up to ^ 20%. 
The 4 simple gravitational lens systems (PGl 115+080, 
SBS1520+530, B1600+434, HE2149-2745) for which time 
delays have been reliably measured yield i/g = 48 ± 3 
km/s/Mpc if the lenses are assumed to have isothermal 
halos of dark matter, while it would be 71 ± 3 km/s/Mpc 
if the lenses instead had constant M / L (Kochanek & 
Schechter (2003). However Koopmans et al. (2003) have 
recently obtained Hq = 75lg km/s/Mpc from a detailed 
reanalysis of the system B1608-I-656, significantly higher 
than their previous estimate of 59 ± 2 km/s/Mpc for this 
system (Fassnacht et al. 2002). We note that Parodi et al. 
(2000) find = 59 ± 4 km/s/Mpc using SN la (see also 
Branch 2000). 

Fig 121 shows a compilation of distances to Type la su- 
pernovae, S-Z clusters and gravitationally lensed systems 
with z > 0.05. The best-fit flat model (not shown) has 
rio = 0.9 ±0.5 and Hq = 60±11 km/s/Mpc. Models with 
= I, Hq = 60 km/s/Mpc and fl^, = 0.3, Qa = 0.7, 
Hq = 72 km/s/Mpc are shown (solid lines); both fit the 
data well. The model with ilm = 1, Hq = 46 km/s/Mpc 
(broken line) is clearly a less good fit to the data. 

In conclusion a Hubble constant in the range 55-65 
km/s/Mpc seems entirely plausible at the present time. 
The value we require, 46 km/s/Mpc, is still below this 
range, but perhaps only by ~ 1 — 2a. We believe the 
present paper provides a powerful stimulus for further 
work on the cosmological distance scale. 

There is however one observational constraint that our 
E-deS model fails to reproduce quite badly, viz. the ampli- 
tude of matter fluctuations on cluster scales {8h~^ Mpc). 
This model has an amplitude of erg 1.1, which is much 
higher than required to match the local abundance of clus- 
ters and weak lensing measurements. The discrepancy is at 
least at the 5 a level, even taking into account the scatter 
in the determinations from different analyses: current es- 
timates for (jg from clusters and from weak lensing range 
from 0.45 to 0.6 (for ilm = 1). This is consistent with 
measurements of peculiar velocity flelds, e.g. the redshift- 
space distortion seen in 2dFGRS on scales < 10h~^ Mpc 
yields /3 = /b = 0.55 ± 0.1 for an E-deS universe, 
where h ^ l/cg is the linear bias parameter (Hawkins 
et al. 2002). * As WMAP cannot probe very small scales, 
one might imagine that further modification of the primor- 
dial spectrum could remove this discrepancy. For example 
Barriga et al. (2001) invoked a 'step' in the spectrum at 
fc ~ 0.1/h Mpc to decrease erg significantly below the cor- 
responding value for a H-Z spectrum, and thus accounted 



* Their best- fit is /3 = 0.49 ±0.09 at the mean redshift of the 
survey: z « 0.17. Combined with the estimate of fe = 1.04 ± 
0.11 from analysis of the bispectrum of 2dFGRS on scales k ~ 
(0.1 - 0.5) /i/Mpc (Verde et al. 2002), this yields D,m(z = 0) = 
0.23±0.09. However, this constraint on the bias, obtained by an 
elaborate statistical analysis, may not be reliable if the biasing 
process is more complicated than envisaged in this work. 
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Fig. 3. CMB angular power spectra for models with 
f^ni = 1, Ho = 46 km/s/Mpc and f^b/i^ 0.019 nor- 
malised to (Tg = 1 with different power law indices for the 
primordial fluctuations. Note the constancy of the ampli- 
tude around £ ~ 900 

for the observed slow evolution of the number density of 
rich clusters with redshift. However, an examination of 
the CMB power spectrum normalized to erg = 1 for various 
power law indices (see Fig.E)) reveals that the implied am- 
plitude on the scale I ~ 900 is essentially constant. Thus 
the recent measurements on this scale by ground-based 
CMB experiments (|Fearson et al. 20021 IKuo et al. 200l 
IRuhl et al. "2001 [Gramge et al. 20031 ), if taken to be reli- 
able, cannot be matched if CTg is made significantly smaller 
than unity by modifying the primordial spectrum. We con- 
clude therefore that an E-deS model with CDM alone can- 
not accommodate both data sets, independently of the 
shape of the primordial spectrum. 

4.2. Modifying the matter content 

In this last section, we examine whether the above dis- 
crepancy can be avoided if we do not restrict ourselves 
to pure CDM models but modify the matter content. 
Specifically, we now have evidence that all the known neu- 
trinos have masses which are rather close to each other, 
with A 171^ ~ 7 X 10~^ eV^ for the electron and muon neu- 
trinos and A ~ 3 x 10""^ eV^ for the muon and tau 
neutrinos, indicated respectively by the oscillation inter- 
pretation of the Solar and atmospheric neutrino anoma- 
lies (see Gonzalez-Garcia & Nir 2003). Moreover the di- 
rect kinematic limit on the neutrino mass from the Mainz 
and Troitsk tritium /3-decay experiments is 2.2 eV (see 
Weinheimer 2002). The addition of massive neutrinos is 
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known to damp the power spectrum on scales smaller than 
their free-streaming length di, — 1230(TOi,/eV)~^ Mpc, 
and thus to lower cg (see Primack & Gross 2000). We 
have therefore introduced 3 quasi-degenerate neutrinos of 
mass 0.8 eV each and find the following model provides 
an acceptable fit: h = 0.46, = ^hh'^ = 0.021, t = 0.10, 
fci = 0.009 Mpc-\ ni = 0.98, = 0.87, fl^ = 0.12. The 
amplitude on clusters scales is perhaps still too high at 
CTg = 0.64, but it is certainly premature to rule out this 
model for this reason alone. 
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Fig. 4. The temperature power spectrum for the best-fit 
power-law ACDM model (dotted black line) from Spergel 
et al. (2003), and for our broken-power-law models (both 
having f^A = 0) with — 0.12 (dot-dashed blue line) 
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In a second model we introduce a small amount of 
'quintessence' with Wq = 0. Such possibility arises natu- 
rally as an attractor solution of a simple exponential po- 
tential: V — Mp exp(— A$/Af — p) which is well motivated 
theoretically (Ratra & Peebles 1988; Wetterich 1988) and 
has been argued to give good agreement with observations 
with ~ 0.1 (Ferreira & Joyce 1997; 1998). Indeed, we 
find that the power at small scales is suppressed in this 
model compared to pure cold dark matter, yielding an ac- 
ceptable fit to the CMB and LSS data (with erg = 0.5) for 
the following parameters: h = 0.45, ojb = ilbh^ — 0.019, 
T = 0.10, ki = 0.012 Mpc-\ 711 = 1.00, 772 = 0.90, 
flQ = 0.12. 

For both cases, the best-fit optical depth to last scat- 
tering T ~ 0.1 is significantly smaller than the value 
of T = 0.17 obtained in the WMAP team's fit to the 



Fig. 5. The temperature-polarization (TE) cross power- 
spectrum and the polarization (EE) power spectrum for 
our E-deS models with f^i^ = 0.12 (dot-dashed blue line) 
and fiq = 0.12 (solid green line), both with r = 0.1, com- 
pared to the concordance ACDM model (dotted line) . The 
thin lines are obtained adopting the higher optical depth 
(r = 0.17) suggested by the WMAP fit to the TE data 
(Kogut et al. 2003). 

temperature-polarization cross-correlation (TE) spectrum 
(Kogut et al. 2003). Nevertheless the predicted TE spec- 
tra for our E-deS models are in reasonable agreement with 
the WMAP data, although the predicted polarization au- 
tocorrelation (EE) spectra differ (see Fig. O . The bary- 
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onic content Wb = i^hh^ — 0.019 is at the upper end of 
the range suggested by considerations of primordial nu- 
cleosynthesis (Fields & Sarkar 2002). The cluster baryon 
fraction expected in these two models is /b ~ 10 — 12% 
(keeping in mind that the 'dark energy' component does 
not cluster on these scales); although low this is consis- 
tent with recent estimates (Sadat & Blanchard 2001). It 
is also necessary to examine the agreement of these mod- 
els with LSS data, in particular the power spectrum ob- 
tained from studies of galaxy clustering and the Lyman-a 
forest. For this purpose we adopt a bias parameter given 
by 6 = 1/(78 ■ As seen in Fig. |H1 both models are then 
in agreement with the APM power spectrum (Baugh & 
Efstathiou 1993; Peacock 1997), the 2dFGRS power spec- 
trum (Percival et al. 2001; Tegmark, Hamilton & Xu 2002) 
and the Lya forest (matter) power spectrum (Croft et 
al. 2002) which has been used recently to set interesting 
constraints on models of large scale structure (Douspis, 
Blanchard & Silk 2001). The model with massive neutri- 
nos provides a particularly good description of LSS data. 
Such a model has already been considered by Elgar0y & 
Lahav (2003) as providing a good fit to the 2dFGRS data, 
but maintaining a constant power law index on large scales 
as they do, then gives a very poor fit to the WMAP data. 
We emphasise that although our model is apparently in 
conflict with the upper bound of il^h^ < 0.0076 quoted 
by Spergel et al. (2003), this latter bound was obtained 
under more restrictive asumptions (in particular adopting 
'priors' on the bias parameter, matter density and Hubble 
parameter), hence is not sufficiently conservative. 

5. Further tests 

It is important to devise means of breaking the Tim — erg 
degeneracy' and establish reliably whether we do live in 
an low matter density universe, and also to devise further 
tests for A independently of the SN la Hubble diagram. 

Attempts to determine J7ni from measurements of pe- 
culiar velocity fields have relied on constraining the bias 
parameter by examination of the statistical properties of 
the galaxy distribution (e.g. Verde et al. 2002). More de- 
tailed examination of the clustering properties of different 
galaxy populations in the ongoing Sloan digital sky survey 
will sharpen this test further (Szapudi et al. 2002). 

Another approach is based on measurement of the 
mean relative peculiar velocity of galaxy pairs as a func- 
tion of their separation (Juskiewicz et al. 2000); a re- 
cent application of this method yields Sim — 0.29to'o9j 
CTg = 0.95lg'^5 without any prior assumptions concern- 
ing the primordial fluctuations or cosmological parameters 
(Feldman et al. 2003). However the inclusion of a compo- 
nent of hot dark matter, as in our model, will affect pecu- 
liar velocities on the relevant small scales, so simulations 
incorporating non-linear effects are necessary to assess the 
robustness of this test (J. Silk, private communication). 

Measurements of 'cosmic shear' induced through grav- 
itational lensing offer yet another possible way to deter- 
mine the matter density and bias separately (see Van 
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Fig. 6. The power spectrum of large scale structure at 
z = and z = 2.3 for our E-deS models with fl^, = 0.12 
(dot-dashed blue lines) and SIq = 0.12 (solid green lines). 
The Lyman-a data have been shifted downwards by 20% 
(corresponding to the la uncertainty in the calibration). 

Waerbeke et al. 2002). Again several ongoing and pro- 
posed large area surveys should allow adequate control of 
systematic uncertainties and enable cosmological param- 
eters to be obtained without prior assumptions. 

A complementary approach is to seek direct evidence 
in the CMB for the presence of a cosmological constant. 
As mentioned earlier, the absence of the expected ISW 
effect at large angular scales in the WMAP data has been 
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disappointing in this respect, but might be regarded as 
a statistical fluctuation. The expected ISW correlations 
with other tracers of large-scale structure are being sought 
but results are not definitive as yet. The power spectrum 
of the E-mode polarization offers an additional way to 
distinguish among models. As seen in Fig.[Sl the nominal 
EE spectra for our E-deS models differ significantly from 
that of the concordance model since the best-fit value of 
r is smaller. However the value of r for our E-deS mod- 
els can be raised to be closer to the WMAP value of 0.17 
(thin lines in Fig. El without significantly affecting the 
Ci or P{k) fits. Hopefully analysis of further data from 
WMAP, as well as other CMB experiments, can distin- 
guish between these possibilities. The most stable differ- 
ence between our E-deS models and the ACDM concor- 
dance model is in fact the matter power spectrum shape in 
the range k ~ (0.01 — 0.03) h/Mpc, which galaxy surveys 
may be able to investigate, provided the possible biasing 
is reliably understood on these scales (Durrer et al. 2003). 

6. Discussion 

We have shown that when the assumption of a single 
power law for the primordial fluctuation spectrum is re- 
laxed, an Einstein-de Sitter model with zero A can fit the 
CMB data as well as if not better than the best ACDM 
concordance model. This is a clear and direct indication 
that the CMB data alone does not require the introduc- 
tion of a non-zero cosmological constant. However a model 
with only cold dark matter cannot simultaneously match 
both the CMB data and the amplitude of matter fluc- 
tuations as indicated by clusters, peculiar velocity flelds 
and weak lensing measurements. We have shown further 
that acceptable Einstein-de Sitter models are indeed pos- 
sible provided they comprise a small amount of a 'dark 
energy' component which does not cluster on small scales, 
such as relic neutrinos with masses of order eV or a pres- 
sureless quintessence field. These models have a low, but 
as we have argued, not unimaginable, Hubble constant. 
Moreover this provides further motivation for laboratory 
experiments sensitive to eV-scale neutrino masses, since a 
detection would provide crucial input for cosmology. 

Given the need to suppress the amplitude of fiuctu- 
ations on cluster scales in any model, we conclude that 
extant CMB and LSS data actually imply the existence of 
a dark component beyond cold dark matter, with a den- 
sity contribution of about 10% of the critical density and 
an equation of state corresponding to conventional pres- 
sureless matter. Such models reproduce quite well the ob- 
served properties of the large scale structure of of universe 
without further adjustment. 

The essential conclusion is that an Einstein-de Sitter 
universe is not yet ruled out, as seems to be generally 
believed. Although there is indeed conflict with some as- 
tronomical observations, we have argued that these data 
are not established beyond reasonable doubt. Given the 
severe coincidence problem associated with a classical cos- 
mological constant as well as the profound implications of 



de Sitter space-time for fundamental physics (e.g. Witten 
2001, Banks & Dine 2001, Dyson, Kleban & Susskind 
2002), it is surely worth investigating these issues further. 
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